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I review constraints on possible New Physics interactions from D — D mixing measurements. 
I consider the most general low energy effective Hamiltonian and include leading order QCD 
running of effective operators. I discuss constraints from an extensive list of popular New 
Physics models, each of which could be discovered at the LHC, that can generate these op- 
erators. In most of the scenarios, strong constraints that surpass those from other search 
techniques could be placed on the allowed parameter space using the existent evidence for 
observation of D meson mixing. 

1 Introduction 

Meson-antimeson mixing has traditionally been of importance because it is sensitive to heavy 
degrees of freedom that propagate in the underlying mixing amplitudes. Estimates of the charm 
quark and top quark mass scales were inferred from the observation of mixing in the K° and B d 
systems, respectively, before these particles were discovered directly. 

This success has motivated attempts to indirectly detect New Physics (NP) signals by com- 
paring the observed meson mixing with predictions of the Standard Model (SM). K°-K° mixing 
has historically placed stringent constraints on the parameter space of theories beyond the SM 
and provides an essential hurdle that must be passed in the construction of models with NP. The 
large mixing signal in the and B s systems, observed at the B-factories and the Tevatron col- 
lider, can be precisely described in terms of the SM alone, which makes the parameter spaces of 
various NP models increasingly constrained. These facts influenced theoretical and experimen- 
tal studies of D° flavor oscillations, where the SM mixing rate is sufficiently small that the NP 
component might be able to compete. There has been a flurry of recent experimental activity 



regarding the detection of D°-D° mixing, which marks the first time Flavor Changing Neutral 
Currents (FCNC) have been observed in the charged +2/3 quark sector. With the potential 
window to discern large NP effects in the charm sector and the anticipated improved accuracy for 
future mixing measurements, the motivation for a comprehensive up-to-date theoretical analysis 
of New Physics contributions to D meson mixing is compelling. 

The phenomenon of meson-anti-meson mixing occurs in the presence of operators that change 
quark flavor by two units Those operators can be generated in both the Standard Model 
and many possible extensions of it. They produce off-diagonal terms in the meson-anti-meson 
mass matrix, so that the basis of flavor eigenstates no longer coincide with the basis of mass 
eigenstates. Those two bases, however, are related by a linear transformation, 

\D 1 )=p\D )± q \D°), (1) 

2 

where the complex parameters p and q are obtained from diagonalizing the D° — D° mass matrix. 
Neglecting CP-violation leads to p = q = l/y/2. The mass and width splittings between those 
mass eigenstates are given by 

X D = — p , VD = • (2) 

It is expected that xq and yo should be rather small in the Standard Model, which is usually 
attributed to the absence of superheavy quarks destroying Glashow-Iliopoulos-Maiani (GIM) 
cancellation. In Eq. ([2]), To is the average width of the two neutral D meson mass eigenstates. 
The quantities which are actually measured in most experimental determinations of the mass 
and width differences, y^ P \ x'^, and y' D , are defined as 



yD° P) = 2/d cos (j)-XB sin - A prod 

x' D = xd cos 5ktt + VD sin 5k-r , (3) 
Ud = Ud cos 5 K n - x D sin 5 Kir , 

where A pro( i = \ N D o — N-^oj / \ Nu° + ^d°) ^ s ^ ne so-called production asymmetry of D° and 

D° (giving the relative weight of D° and D° in the sample) and Sk-k is the strong phase difference 
between the Cabibbo favored and double Cabibbo suppressed amplitudes^!, which is usually 
measured in D — > Kir transitions. In what follows we shall negle ct CP-violating parameters (j) 
and A m . In this limit = Ud- Please see recent reviews ^ - ^^ for more complete analysis. 

2 Experimental Constraints on Charm Mixing 

The recent interest in D°-D° mixing started with the almost simultaneous observations by the 
BaBarEl and Belled collaborations of nonzero mixing signals at about the per cent level, 

y' D = (0.97 ± 0.44 ± 0.31) • 10~ 2 (BaBar) , (4) 
yf P) = (1.31 ± 0.32 ± 0.25) • 10~ 2 (Belle) . (5) 

This was soon followed by the announcement by the Belle collaboration of mixing measurements 
from the Dalitz plot analyses of D° — > iTs7r + 7r - ^, 

XD = (0.80 ±0.29 ±0.17) • 10~ 2 , y D = (0.33 ± 0.24 ± 0.15) • 10" 2 . (6) 

A fit to the current database by the Heavy Flavor Averaging Group (HFAG) gives^ 

x D = 9.8t 2 : ? • 10" 3 , y D = (7.5 ±1.8)- 10~ 3 , (7) 
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which is obtained assuming no CP-violation affecting mixing. It is important to note that the 
combined analysis of xd and ?/d excludes the "no-mixing" point xd = y~£> = by 6.7a ^ . This 
fact adds confidence that charm mixing has indeed been observed. Then, a correct interpretation 
of the results is important. In addition, as with any rare low-energy transition, the question 
arises on how to use it to probe for physics beyond the Standard Model. 

3 Standard Model "background" in D° — D° mixing 

Theoretical predictions for xd and yr> obtained in the framework of the Standard Model histor- 
ically span several orders of magnitude. I will not discuss predictions of th e SM for the charm 
mixing rates here, instead referring the interested reader to recent reviews HEBl it might be 
advantageous to note that there are two approaches to describe D° — D° mixing, neither of 
which give very reliable results because m c is in some sense intermediate between heavy and 

light. r-r-i 

The inclusive approach USUI j s based on the operator product expansion (OPE). In the 
formal limit m c ^> A limit, where A is a scale characteristic of the strong interactions, xd and 
yo can be expanded in terms of matrix elements of local operators. The use of the OPE relies 
on local quark- hadron duality, and on A/m c being small enough to allow a truncation of the 
series after the first few terms. This, however, is not realized in charm mixing, as the leading 
term in l/m c is suppressed by four and six powers of the strange quark mass for xd and yp 
respectively. The parametrically-suppressed higher order terms in l/m c can have less powers 
of m s , thus being more important numerically^! This results in reshuffling of the OPE series, 
making it a triple expansion in l/m c , m s , and a s . The (numerically) leading term contains over 
twenty matrix elements of dimension-12, eight-quark operators, which are difficult to compute 
reliably. A naive power counting then yields xd,Vd < 10 -3 . The exclusive approach ^21 sums 
over intermediate hadronic states. Since there are cancellations between states within a given 
SU(3) multiplet, one needs to know the contribution of each state with high precision. However, 
the D is not light enough that its decays are dominated by a few final states. In the absence 
of sufficiently precise data, one is forced to use some assumptions. Large effects in yp appear 
for decays close to D threshold, where an analytic expansion in SU(3)f violation is no longer 
possible. Thus, even though theoretical calculations of x n a nd yn are quite uncertain, the values 
^d~2/d~1% are quite natural in the Standard Model^. 

It then appears that experimental results of Eq. ([7]) are consistent with the SM predic- 
tions. Yet, those predictions are quite uncertain to be subtracted from the experimental data to 
precisely constrain possible NP contributions. In this situation the following approach can be 
taken. One can neglect the SM contribution altogether and assume that NP saturates the result 
reported by experimental collaborations. This way, however, only an upper bound on the NP 
parameters can be placed. A subtlety of this method of constraining the NP component of the 
mixing amplitude is related to the fact that the SM and NP contributions can have either the 
same or opposite signs. While the sign of the SM contribution cannot be calculated reliably due 
to hadronic uncertainties, xd computed entirely within a given NP model can be determined 
rather precisely. This stems from the fact that NP contributions are generated by heavy degrees 
of freedom making short-distance OPE reliable. This means that only the part of parameter 
space of NP models that generate X£> of the same sign as observed experimentally can be reliably 
and unambiguously constrained. 

4 New Physics contributions to D° — D° mixing 

Any NP degree of freedom will generally be associated with a generic heavy mass scale M, 
at which the NP interaction will be most naturally described. At the scale m c of the charm 



mass, this description will have been modified by the effects of QCD, which should be taken into 
account. In order to see how NP might affect the mixing amplitude, it is instructive to consider 
off-diagonal terms in the neutral D mass matrix, 
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where the first term contains TL^ C= ~ 2 , which is an effective |AC| = 2 hamiltonian, represented 
by a set of operators that are local at the ~ mjj scale. Note that a 6-quark also gives a 
(negligible) contribution to this term. This term only affects xd, but not yrj. 

The second term in Eq. ([5]) is given by a double insertion of the effective |AC| = 1 Hamil- 
tonian TL^ C= ~ 1 . This term is believed to give dominant contribution to D° — D° mixing in the 
Standard Model, affecting both x and y. It is generally believed that NP cannot give any siz- 
able contribution to this term, since Ti!^ c= ~ 1 Hamiltonian also mediates non-leptonic D-decays, 
which should then also be affected by this NP contribution. I will show that there is a well- 
defined theoretical limit where NP contribution dominates lifetime difference yrj and consider 
implications of this limit in " real world" . 

4-1 New Phyiscs in |AC| = 1 interactions. 

Consider a non-leptonic D° decay amplitude, A[D° — > n], which includes a small NP contri- 
bution, A[Z) — > n] = An M ^ + A^ P \ Here, A^ F ^ is assumed to be smaller than the current 
experimental uncertainties on those decay rates. This ensures that NP effects cannot be seen in 
the current experimental analyses of non-leptonic D-decays. One can then write yj) as 

vd - Er L 4 SM) 4 SM) + 2Ef L ^ NP) ^ SM) • ( 9 ) 

n D n D 

The first term of Eq. (schematic) represents the SM contribution to t/d- The SM contribution to 
yD is known to vanish in the limit of exact flavor SU(3). Moreover, the first order correction is 
also absent, so the SM contribution arises only as a second order effect^! This means that in the 
flavor SU(3) limit the lifetime difference yv is dominated by the second term in Eq. ([9]), i.e. New 
Physics contributions, even if their contibutions are tiny in the individual decay amplitudes^! 
A calculation reveals that NP contribution to y^ can be as large as several percent in R-parity- 
violating SUSY mo delsEl or as small as ~ 10 10 in the models with interactions mediated by 
charged Higgs particles^. 

This wide range of theoretical predictions can be explained by two observations. First, many 
NP affecting |AC| = 1 transitions also affect \AB\ = 1 or | A.S*| = 1 decays or kaon and B-meson 
mixings, which are tightly constrained. Second, a detailed look at a given NP model that can 
potentially affect yo reveals that the NP contribution itself can vanish in the flavor SU(3) limit. 
For instance, the structure of the NP interaction might simply mimic the one of the SM. Effects 
like that can occur in some models with extra space dimensions. Also, the chiral structure 
of a low-energy effective lagrangian in a particular NP model could be such that the leading, 
mass-independent contribution vanishes exactly, as in a left-right model (LRM). Finally, the NP 
coupling might explicitly depend on the quark mass, as in a model with multiple Higgs doublets. 
However, most of these models feature second order S*C/(3)-breaking already at leading order 
in the l/m c expansion. This should be contrasted with the SM, where the leading order is 
suppressed by six powers of m s and term of order m 2 s only appear as a 1/m^-order correction. 

4-2 New Phyiscs in |AC| = 2 interactions. 



Though the particles present in models with New Physics may not be produced in charm quark 
decays, their effects can nonetheless be seen in the form of effective operators generated by the 



exchanges of these new particles. Even without specifying the form of these new interactions, 
we know that their effect is to introduce several |AC| = 2 effective operators built out of the 
SM degrees of freedom. 

By integrating out new degrees of freedom associated with new interactions at a scale M, 
we are left with an effective hamiltonian written in the form of a series of operators of increasing 
dimension. Operator power counting then tells us the most important contributions are given 
by the operators of the lowest possible dimension, d = 6 in this case. This means that they 
must contain only quark degrees of freedom and no derivatives. Realizing this, we can write 
the complete basis of these effective operators, which can be done most conveniently in terms 
of chiral quark fields, 

(f\n NP \i) = G^difi) (f\Qi\i)(n) , (10) 

i=l 

where the prefactor G has the dimension of inverse-squared mass, the Cj are dimensionless 
Wilson coefficients, and the Qi are the effective operators: 



Q2 = (ul^cl) (ur^cr) , Q G = (ur-j^cr) {ur^cr) , 

Q3 = (u L c R ) (urc l ) , Q 7 = (u l cr) {u l c r ) , 

Qi = {urc l ) {urc l ) , Qs = {ul^^cr) {ul^cr) 



(11) 



In total, there are eight possible operator structures that exhaust the list of possible independent 
contributions to |AC| = 2 transitions. Since these operators are generated at the scale M where 
the New Physics is integrated out, a non-trivial operator mixing can occur when one takes 
into account renormalization group running of these operators between the scales M and fi, 
with fi being the scale where the hadronic matrix elements are computed. We shall work at 
the renormalization scale [i = m c ~ 1.3 GeV. This evolution is determined by solving the RG 
equations obeyed by the Wilson coefficients, 

" C(^)=7 T C(/i) , (12) 



(flog /i 



where 7 represents the matrix of anomalous dimensions of the operators in Eq. (jlipH^l. Due to 
the relatively simple structure of 7, one can easily write the evolution of each Wilson coefficient 
in Eq. (jlOf) from the New Physics scale M down to the hadronic scale u, taking into account 
quark thresholds. Corresponding to each of the eight operators {Qi} (i = 1, . . . , 8) is an RG 
factor ri{fi,M). The first of these, ri(/i,M), is given explicitly by 



\a s (m t )J \a s (m b )J \ a s (u) 



6/25 



(13) 



and the rest can be expressed in terms of 7*1 {fj,, M) as 
r 2 {u,M) = [r l { f i,M)} 1 / 2 , 

r 3 {u, M) = [ n {fi, M)]- 4 , r 6 {fi, M) = n(/i, M) 

r 4 {a, M) = [n(M, M)]( 1 +v / 24T)/6 ; ^ M ) = r ^ M ) 

r 5 {u, M) = [n(/i, M)]^-^ 21 )/ 6 , r s (p, M) = r 5 (/i, M) 



(14) 



The RG factors are generally only weakly dependent on the NP scale M since it is taken to be 
larger than the top quark mass, nit, and the evolution of a s is slow at these high mass scales. In 
Table[H we display numerical values for the ri{fi, M) with M = 1, 2 TeV and \x = m c ~ 1.3 GeV. 
Here, we compute a s using the one-loop evolution and matching expressions for perturbative 
consistency with the RG evolution of the effective hamiltonian. A contribution to D° — D° 



M(TeV) 


ri (m c ,M) 


r 2 (m c ,M) 


r 3 (m c ,M) 


r 4 (m c ,M) 


r 5 (m c ,M) 


1 


0.72 


0.85 


3.7 


0.41 


2.2 


2 


0.71 


0.84 


4.0 


0.39 


2.3 



Table 1: Dependence of the RG factors on the heavy mass scale M. 



mixing from a particular NP model can be obtained by calculating matching conditions for the 
Wilson coefficients Cj at the scale M, running their values down to \x and computing the relevant 
matrix elements of four-quark operators. A generic model of New Physics would then give the 
following contribution x^, 



G 



fl)B D m D 



2 5 7 

~[Ci(m c ) + C 6 (m c )] - — [C 4 (m c ) + C 7 {m c ]) + —C 3 (m c 



5C 2 (m c ) 
6 



+ [C 5 (m c ) + C 8 (m c )\ 



(15) 



Here we simplified the result by assuming that all non-perturbative ('bag') parameters are equal 
to .Bd ~ 0.82. The Wilson coefficients at the scale \x are related to the Wilson coefficients at 
the scale M by renormalization group evolution, 



Ci(m c 
C 2 (m c 

C 3 (m c 
Carrie 



C 5 (m c ) 

C 6 (m c ) 
C 7 (m c ) 

C 8 (m c ) 



n(m c ,M)Ci(M) , 
r 2 (m c ,M)C 2 (M) , 

I [r 2 (m c , M) - r 3 (m c , M)\ C 2 {M) + r 3 (m c , M)C 3 (M) 



x 241 



[r 5 (m c , M) - r 4 (m c , M)] 
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C 4 (M) + -C b {M) 



+ - [r 4 (m c ,M)+r 5 (m c ,M)} C 4 (M) 



1 



8V241 



[r 4 (m c , M) - r 5 (m c , M)\ [C 4 (M) + 64C 5 (M)] 



1 



+ - [r 4 (m c ,M)+r 5 (m c ,M)] C 5 (M) 



r 6 (m c ,M)C 6 (M) , 

[r 8 (m c ,M)-r 7 (m c ,M)} 



'2<il 
1 



(16) 
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C 7 (M) + -Cs(M) 



+ - [r 7 (m c , M) + r 8 (m c , M)] C 7 (M) 



1 



8V241 



[r 7 (m c , M) - r 8 (m c , M)] [C 7 (M) + 64C 8 (M)] 



1 



+ - [r 7 (m c , M) + r 8 (m c , M)] C 8 (M) , 

A contribution of each particular NP model can then be studied using Eq. (|15p . Even before 
performing such an analysis, one can get some idea what energy scales can be probed by D° — D° 
mixing. Setting G = 1/M 2 and Cj(M) = 1, we obtain M ~ 10 3 TeV. More realistic models 
can be probed in the region of several TeV, which is very relevant for LHC phenomenology 
applications. 

A program described above has been recently executed^! for 21 well-motivated NP models, 
which will be actively studied at LHC. The results are presented in Table [2 As can be seen, out 
of 21 models considered, only four received no useful constraints from D° — D° mixing. More 
informative exclusion plots can be found in that paper^^l as well. It is interesting to note that 



Model 


Approximate Constraint 


Fourth Generation 


\V,, h iVrh>\ ■ my < 0.5 (GeV) 

| ' (10 * CO | ^ " * " \ w y y 


Q = —1/3 Singlet Quark 


s 2 ■ m s < 0.27 (GeV) 


Q = +2/3 Singlet Quark 


|A„J < 2.4 • 10~ 4 


Little Higgs 


Tree: See entry for Q = —1/3 Singlet Quark 




Box: Parameter space can reach observed xd 


Generic Z' 


Mz'/C > 2.2 • 10 3 TeV 


Family Symmetries 


mi// > 1.2 • 10 3 TeV (with mi/m 2 = 0.5) 


Left-Right Symmetric 


No constraint 


Alternate Left-Right Symmetric 


Mr > 1.2 TeV (m Dl = 0.5 TeV) 




(Am/mn-,)/M R > 0.4 TeV" 1 


Vector Leptoquark Bosons 


Mvto > 55(A P p/0.1) TeV 


Flavor Conserving Xwo-Higgs-Doublet 


No constraint 


Flavor Changing Neutral Higgs 


m H /C > 2.4 • 10 3 TeV 


FC Neutral Higgs (Cheng-Sher) 


77iH-/|A uc > 600 GeV 


Scalar Leptoquark Bosons 


See entry for RPV SUSY 


Higgsless 


M > 100 TeV 


Universal Extra Dimensions 


No constraint 


Split Fermion 


M/\Ay\ > (6 • 10 2 GeV) 


Warped Geometries 


Mi > 3.5 TeV 


MSSM 


|(Si 2 )lr,Rl| < 3.5 • 10~ 2 for rh ~ 1 TeV 




I(5i2)ll,rr| <-25 form- 1 TeV 


SUSY Alignment 


m > 2 TeV 


Supersymmetry with RPV 


A'mA'nfc/m^ < 1.8 • 10- 3 /100 GeV 


Split Supersymmetry 


No constraint 



Table 2: Approximate constraints on NP models from D° mixing. 



some models require large signals in the charm system if mixing and FCNCs in the st range and 
beauty systems are to be small (as in, for example, the SUSY alignment 



5 Conclusions 



I reviewed implications of recent measurement of D° — D mixing rates for constraining models 
of New Physics. A majority of considered models received competitive constraints from D° — D° 
mixing measurements despite hadronic uncertainties that plague SM contributions. It should be 
noted that vast majority of predictions of NP models do not suffer from this uncertainty, and 
can be computed reliably, if lattice QCD community provides calculations of matrix elements of 
four-fermion operators Eq. (|11|) . 

Another possible manifestation of new p hysics i nteractions in the charm system is associated 
with the observation of (large) CP-violation 1H4H3 This is due to the fact that all quarks that 
build up the hadronic states in weak decays of charm mesons belong to the first two generations. 
Since 2x2 Cabbibo quark mixing matrix is real, no CP-violation is possible in the dominant tree- 
level diagrams which describe the decay amplitudes. CP-violating amplitudes can be introduced 
in the Standard Model by including penguin or box operators induced by virtual 6-quarks. 
However, their contributions are strongly suppressed by the small combination of CKM matrix 
elements V^V^,. It is thus widely believed that the observation of (large) CP violation in charm 
decays or mixing would be an unambiguous sign for New Physics. 
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